Panel Stiffener Debonding Analysis Using a Shell/3D Modeling Technique 


Ronald Krueger *, James G. Ratcliff e*, Pierre J. Minguet** 
* National Institute of Aerospace, **The Boeing Company 


Abstract 

A shear loaded, stringer reinforced composite panel is analyzed to evaluate the fidelity of 
computational fracture mechanics analyses of complex structures. Shear loading causes the 
panel to buckle. The resulting out-of-plane deformations initiate skin/stringer separation at 
the location of an embedded defect. The panel and surrounding load fixture were modeled 
with shell elements. A small section of the stringer foot, web and noodle as well as the 
panel skin near the delamination front were modeled with a local 3D solid model. Across 
the width of the stringer foot, the mixed-mode strain energy release rates were calculated 
using the virtual crack closure technique. A failure index was calculated by correlating the 
results with a mixed-mode failure criterion of the graphite/epoxy material. The objective 
was to study the effect of the fidelity of the local 3D finite element model on the computed 
mixed-mode strain energy release rates and the failure index. 

Keywords: A: Structural composites; B: Delamination, C: Finite element analysis (FEA), 
C: Buckling; Computational fracture mechanics. 


1 Introduction 

Many composite components in aerospace structures are made of flat or curved panels with 
co-cured or adhesively bonded frames and stiffeners. Recent studies focused on the 
investigation of the debonding mechanism and included testing of skin/stiffener panels and 
failure analysis using shell models [1,2]. Over the last decade, a consistent step-wise 
approach has been developed which uses experiments to determine the failure mechanism, 



computational stress analysis to determine the location of first matrix cracking and 
computational fracture mechanics to investigate the potential for delamination growth. 
Testing of thin-skin stiffened panels designed for use in pressurized aircraft fuselages has 
shown that bond failure at the tip of the frame flange is an important and very likely failure 
mode [3]. Comparatively simple, laboratory size coupon type specimens consisting of a 
stringer flange bonded onto a skin have been developed to study the fundamental 
mechanisms of skin/stiffener debonding [4], The failure that initiates at the tip of the flange 
in these coupon type specimens is nearly identical to the failure observed in the full-scale 
panels and frame pull-off specimens [4-6] . A methodology based on fracture mechanics [7] 
has been used successfully to investigate the onset and growth of delaminations in simple 
characterization specimens and laboratory size coupon type specimens [5,6]. Future 
acceptance of the methodology by industry and certification authorities, however, requires 
the successful demonstration of the methodology on a structural level . 

For the demonstration of the methodology on the structural level, a stringer stiffened panel, 
as shown in Fig. 1, had been analyzed previously [8, 9]. The square (1016 mm x 1016 mm) 
panel made of carbon/epoxy tape is reinforced with three stringers made of carbon/epoxy 
plain weave fabric. Details will be discussed in the following sections. During 
manufacturing, an artificial defect of about 82 mm in length was placed at the termination 
of the center stiffener. The stiffened panel was bolted to a steel picture frame and subjected 
to shear loading which caused the panel to buckle as shown in the finite element model 
depicted in Fig. 1 [8,9]. The resulting out-of-plane deformation caused skin/stringer 
separation to initiate at the location of the artificial defect. A small section of the stringer 
foot and the panel skin in the vicinity of the embedded defect were modeled with a local 
3D solid model as shown in the enlargement in Fig. 1. The mixed-mode strain energy 
release rates were calculated using the virtual crack closure technique [10, 11] across the 
width of the stringer foot. A failure index was calculated by correlating the results with the 
mixed-mode failure criterion of the graphite/epoxy material [8,9]. 

The objective of the current research was to study the effect of the fidelity of the local 3D 
finite element model on the computed mixed-mode strain energy release rates and the 



failure index. In the original model, shell elements represented the skin and stiffener foot 
[8,9]. In previous studies, a model was introduced where only the stiffener foot and skin 
were included in a local 3D model, while the web and hat were modeled with shell 
elements as shown in Fig. 1 . The pair of 2D illustrations in Fig. 2 are cross sections of the 
local 3D finite element model [8,9]. For the current study, a new set of models was 
introduced, as shown in Fig. 3, where the stringer web was included in the local solid 
model. Another set of models included the transition radius between the web and foot as 
well as the detailed noodle region in the local solid model as shown in Fig. 4. Details of the 
models are discussed later. Thus, this study considers the effects of various geometrical 
simplifications that might be considered by an analyst when modeling debonding in 
complex structures. The work discussed herein is intended to be used as a guide to 
determine the trade off between modeling complexity and fidelity on the corresponding 
results. The current research complements previous studies [8,9]. 


2 Methodology 

2.1 Linear Elastic Fracture Mechanics 

Linear elastic fracture mechanics has proven useful for characterizing the onset and growth 
of delamination in composite laminates [7]. When using fracture mechanics, the total strain 
energy release rate, G r , is calculated along the delamination front. The term, G T , consists of 
three individual components. The first component, Gj, arises due to interlaminar tension. 
The second component, G n , arises due to interlaminar sliding shear (shear stresses parallel 
to the plane of delamination and perpendicular to the delamination front). The third 
component, G in , arises due to interlaminar scissoring shear (shear stresses parallel to the 
plane of delamination and parallel to the delamination front). The calculated G„ G,„ and 
G m components are then compared to interlaminar fracture toughness values in order to 
predict delamination onset or growth. The interlaminar fracture toughness values are 
determined experimentally over a range of mode mixity from pure mode I loading to pure 
mode II loading [12-14], 



A quasi static mixed-mode fracture criterion is determined by plotting the interlaminar 
fracture toughness, G c , versus the mixed-mode ratio, G„l G T . The fracture criteria is 
generated experimentally using pure Mode I (G„l G T =0) Double Cantilever Beam (DCB) 
tests, pure Mode II (G„/ G T = 1) four point End-Notched Flexure (4ENF) tests, and Mixed- 
Mode Bending (MMB) tests of varying ratios of G, and G„. A fracture criterion was 
suggested by Benzeggah and Kenane [15] using a simple mathematical relationship 
between G c and G„! G T 



as shown in Fig. 5. 


In this expression, G, c and G„ c are the experimentally-determined fracture toughness data 
for mode I and II. The factor rj is determined by a curve fit. Fracture initiation is expected 
when, for a given mixed-mode ratio G„l G T , the calculated total energy release rate, G r , 
exceeds the interlaminar fracture toughness, G c . 

Although several specimens, including the edge-cracked torsion specimen (ECT), have 
been suggested for the measurement of the mode III interlaminar fracture toughness 
property [16, 17], an interaction criterion incorporating the mode III shear, however, has 
not yet been established. Therefore, a modified definition is introduced for three- 
dimensional analysis, which also yields results for the scissoring mode G m . In the modified 
definition, G s denotes the sum of the in-plane shearing components G„ + G,„ [8,9]. This is 
necessary since a mixed-mode failure criterion, which accounts for all three modes, is 
currently not available. For analyses where G/// =0 , this definition is equal to the commonly 
used definition of the mixed-mode ratio, G,J G T mentioned above. 

To determine failure along the delamination front, the critical energy release rate, G c , is 
calculated using Eq. (1) with G„ = G s at each point along the delamination front. 
Subsequently, the failure index GjJ G c is determined from the computed total energy 
release rate, G r , and the critical energy release rate, G c , with the assumption that 
delamination propagation occurs for 




22 Analysis Tools 

22.1 Virtual Crack Closure Technique 

A variety of methods are used in the literature to compute the strain energy release rate 
based on results obtained from finite element analysis. For delaminations in laminated 
composite materials where the failure criterion is highly dependent on the mixed-mode 
ratio, the virtual crack closure technique (VCCT) [10,11] has been most widely used for 
computing energy release rates. VCCT calculations using continuum (2D) and solid (3D) 
finite element analyses provide the mode separation required when using the mixed-mode 
fracture criterion. 

222 A Global/Local Shell 3D Modeling Technique 

Built-up structures are traditionally modeled and analyzed using plate or shell finite 
elements, as shown in Fig. 1, to keep the modeling and computational effort affordable. 
The fidelity of the resulting computed mixed-mode strain energy release rate components, 
however, depends on many variables. These variables include element type, order of the 
shape functions and shear deformation assumptions, kinematic constraints in the 
neighborhood of the delamination front, and continuity of material properties and section 
stiffnesses in the vicinity of the debond when delaminations or debonds are modeled with 
plate or shell finite elements [18]. These problems may be avoided by using three- 
dimensional models. However, since many layers of brick elements through the thickness 
are often necessary to model the individual plies, the size of finite element models required 
for accurate analyses may become prohibitively large. 

For detailed modeling and analysis of the delaminations, the shell/3D modeling technique 
will reduce the modeling time compared to that required to run a fully three-dimensional 
finite element model. The technique will also reduce computational time because only a 
relatively small section of the mesh needs to be modeled with solid elements, minimizing 



the overall size of the model. The technique combines the accuracy of the full three- 
dimensional solution with the computational efficiency of a plate or shell finite element 
model. The technique has been demonstrated for various applications such as fracture 
toughness characterization specimens [19], on the coupon level for the skin/stringer 
separation specimen [20] and in related studies for skin/stringer separation [8, 9]. The 
enlargement in Fig. 1 illustrates the regions within the stringer stiffened panel that are 
modeled with shells and solid elements. 


3 Finite Element Modeling 

In the current study, a finite element analysis of the three-stringer panel shown in Fig. 1 
was conducted. The load frame and the three-stringer panel were modeled with beam and 
shell elements. A small section of the stringer and the panel skin in the vicinity of the 
embedded defect was modeled with a local 3D model. For all analyses, the stiffener hat 
was modeled with shell elements. For all analyses, the finite element software ABAQUS® 
was used which is manufactured by Dassault Systemes Simulia Corp. (DSS), Providence, 
RI, USA. 

3.1 Global Shell Model of Stringer Stiffened Panel 

The global model included the steel load frame and attachments, the panel made of 
graphite/epoxy prepreg tape, and the stringers made of graphite/epoxy fabric, as shown in 
Fig. 1. The outer steel load frame and the attachment bolts were modeled with beam 
elements (ABAQUS® element type B21) [21], The inner steel load frame, which overlaps 
the panel edge, was modeled with shell elements (ABAQUS® element type S4). The 
stiffener components, such as the foot, web and hat were also modeled with shell elements. 

The panel skin and the stiffener foot are modeled as separate entities. The S4 shell 
elements are located at the panel skin and stiffener foot respective mid-planes. The shell 
elements are connected by beam elements (ABAQUS® element type B31 modeled as steel) 
to provide point-to-point constraints between the two surfaces [22] . In the section 
containing the artificial defect, the beam elements were replaced by gap elements 



(ABAQUS® element type GAPUNI). The gap elements allow the modeling of contact 
between two nodes. The nodes can be in contact (gap closed), which prevents element 
interpenetration or separated (gap open) which allows the skin/stringer separation [21]. 

3 2 Local 3D Model of the Stringer Foot and Panel Skin 

The local 3D model of the stringer foot and panel skin was generated using brick elements 
(ABAQUS® element type C3D8I) and consisted of an intact section and a delaminated 
section with a fine mesh around the delamination front. Two examples with different mesh 
refinement (designated BF-8/10/8 and BF-20/20/20) are shown in Fig. 2. The artificial 
defect is located at the bondline between the stringer foot and the panel as shown in the 
enlargements of Fig. 2. This defect was treated as a delamination and modeled as a discrete 
discontinuity using two unconnected nodes with identical coordinates on each side of the 
delamination. Contact was modeled between the delaminated surfaces to avoid 
interpenetration during the analysis. 

Four elements through the thickness were used to model the foot of the stiffener made of 
carbon/epoxy fabric as shown in the 2D cross sections in Fig. 2. The skin ply made from 
carbon/epoxy tape which is adjacent to the plane of the delamination was modeled with 
one element through the thickness. The remaining 10 plies of carbon/epoxy tape were 
modeled with three elements through the thickness as shown in Fig. 2. 

33 Local 3D Model of the Panel Skin, Stringer Foot and Web 

The modeling fidelity was increased by including the stiffener web in the local 3D model 
as shown in Fig. 3 (TFL3-8/10/8, TFL3-20/20/20). The local 3D model shown consisted of 
two separately meshed sections: the stiffener web and the panel skin/stiffener foot. The 
web was modeled with eight brick elements (ABAQUS® element type C3D8I) through the 
thickness. The two models are connected to each other using multipoint constraints 
(ABAQUS® MPC option *TIE). The multipoint constraints are used to tie two surfaces 
together for the duration of a simulation, which makes the translational and rotational 
motion equal for a pair of surfaces. Nodes on the slave surface (defined by the user) are 



constrained to have the same motion as the point on the master surface (defined by the 
user) to which it is closest [21]. 

3.4 Local 3D Model of the Panel Skin, Stringer Foot and Web Including the Noodle 
and Transition Radius 

The modeling fidelity was increased further by including the noodle region and transition 
radius in the local 3D model as shown in Fig. 4 (TN1L3C-8/10/8, TN1L3C-20/20/20). The 
triangular region underneath the web, where the stiffener connects to the skin, is usually 
referred to as a noodle. As above, the local 3D model shown consisted of two separately 
meshed sections which are connected with tie constraints: The model of the stiffener 
including the web, noodle region and transition radius and the model of panel skin and 
delaminated interface is shown in the cross section of Fig. 4. The two solid models were 
joined using the *TIE option in ABAQUS® [21]. The solid model of the stringer foot and 
delamination consisted of brick elements (ABAQUS® element type C3D8I) and included 
an intact section and a delaminated section with a fine mesh around the delamination front, 
as discussed above. 

Solid models representing the stiffener foot, the web and the noodle region consisted of 
eight-node C3D8I brick elements, with a small number of six-node C3D6 prism elements 
used to model part of the noodle region as shown in the cross sections of Fig. 4. Since the 
exact transition radius was not known, the influence of detailed local 3D modeling on 
computed strain energy release rates was studied. In this case, a number of models were 
generated with different foot/web transition radii. Four models were generated with radius 
values of 0.254 mm, and 0.71 1 mm, 2.54 mm and 5.08 mm. Models with radii 2.54 mm 
and 0.71 1 mm were thought to correspond to values used in manufactured panels. The 
small radius 0.254 mm (cross section shown in Fig. 4c) was selected to determine if the 
computed results were similar to those obtained from the model discussed in section 3.3 
which did not included the radius and noodle. The larger radius 5.08 mm (cross section 
shown in Fig. 4c) was chosen as an upper limit. It did not appear to be meaningful from a 
design standpoint to assume larger radii. 



3.5 Finite Element Model Assembly, Load and Boundary Conditions 

The local 3D models were placed to complete the assembled model as shown in the 
enlargement of Fig. 1 . At the boundaries, the nodes along the shell edges were used to 
connect the shell model with the local 3D solid models. The connection was accomplished 
using the shell to solid coupling option in ABAQUS®, which allows the connection 
between non-conforming shell and solid models. The coupling option uses a set of 
internally defined distributing coupling constraints to couple the motion of a row of nodes 
along the edge of the shell model to the motion of a set of nodes defined on a surface of the 
solid model [21]. 

For modeling the experiment, which was performed under displacement control, uniform 
displacements u, v were applied at one comer node to introduce shear as shown in Fig. 1 . 
The in-plane displacements, u, v, were suppressed at the diagonally opposite comer and the 
out-of-plane displacements, w, were suppressed along all four edges across the entire width 
of the inner and outer steel load frame. 

3.6 Analysis Overview 

A total of 41 different model combinations were analyzed. An overview of all models 
generated is given in reference [23]. The following cases are discussed here: 

• Models were analyzed where only the stiffener foot and skin were included in a 
local 3D model, while the web and hat were modeled with shell elements as shown 
in Fig. 2. 

• A set of models were created, as shown in Fig. 3, where the stringer web was 
included in the local solid model. 

• Another set of models included the transition radius between the web and foot as 
well as the detailed noodle region in the local solid model as shown in Fig. 4. The 
stiffener hat, modeled with shell elements, was kept unchanged. Four different 
transition radii were studied (r=0.254 mm, 0.711 mm, 2.54 mm, 5.08 mm). 

• In order to capture the local failure near the edges, models with a locally refined 
fine mesh were chosen as shown in the cross sections of Fig. 2 to 4. In order to 



capture the local failure index distribution in the vicinity of the web termination, 
models with a fine mesh in the center were chosen (also shown in Fig. 2 to 4). 


4 Analysis Results 

4.1 Model Deformation 

Under the applied shear loading, the analysis predicts the buckling deformation shown in 
Fig. 1 . For the simulated delamination length (a=81 .9 mm), three peaks and one trough can 
be observed in the panel bays adjacent to the center stiffener as shown in Fig. 1 . More 
details may be found in a previous study [8] . 

Details of the deformed finite element models are shown in the enlargements of Fig. 2 to 4 
after the entire external displacement w=v=6.35 mm had been applied. Mode I opening was 
observed across the entire width of the stringer over the entire delaminated length. 

4.2 Comparison of Failure Indices Computed from Different Local 3D Models 

For each nodal point along the delamination front, the critical energy release rate, G c , was 
calculated from a mixed-mode failure criterion (Eq. 1) for the computed mixed-mode ratio, 
G s / G T . Subsequently, the failure index G T / G c was calculated from Eq. 2. The failure 
index was calculated for the final load increment and plotted versus the location, s, across 
the width of the stringer, b, as shown in Fig. 1 and 2. 

The computed failure index distributions across the width of the stringer obtained from 
different models were plotted in Fig. 6 to 12 for comparison with reference results. The 
combined results from the models with refined edges and center (Fig. 2) were used as 
reference results in all the figures and are depicted as blue circles. For all cases shown, the 
failure index peaked at the edges (£=0.0 and s=l .0) with an additional peak around the 
center (£=0.5) underneath the stringer web. 


First, the results obtained from a model that included the panel skin, stiffener foot and the 
web in the local solid model as shown in Fig. 3 were plotted in Fig. 6 for comparison with 



the reference result. Qualitatively, both distributions followed the same trend. In two areas 
to the left and right of the web (0.0 £ 0.4 and 0.7 s s s 1 .0), the results are in good 

agreement. Locally, near the web (0.4 s ss 0.7), the distributions differ. The peak values 
computed for the local solid model that included the web (red squares) exceed the 
reference values (blue circles) by about 49%. The peak location is also offset. 

The results obtained for a transition radius r=0.71 1 mm are plotted in Fig. 7. In the area to 
the left of the web (0.0 £ sz 0.4), the results are in good agreement. In the areas to the right 
of the web (0.7 s s s 1 .0), the results are higher for the models which included the web and 
the noodle (green diamonds). An additional small plateau was observed for 0.47 < s < 0.49 
(G r / G c =29) which was not observed in the results used as reference. Also, the peak 
values computed for the local solid model that included the web and the noodle exceed the 
reference values (blue circles) by about 67%. As before, the peak location is offset 
compared to the reference results. 

For the results computed for the larger transition radius (r=2.54 mm), an additional local 
maximum was observed for 0.44 < s < 0.46 ( G T / G c ~ 15) as shown in Fig. 8. In the areas 
to the right of the web (0.7 £$£ 1.0), the results are higher for the models, which included 
the web and the noodle (orange triangles). Also, the peak values computed for the local 
solid model that included the web and the noodle exceed the reference values (blue circles) 
by about 51%. For this transition radius, the peak location is almost identical to the peak 
observed for the reference results. 

The results obtained for a smaller transition radius (r=0.254 mm as shown in Fig. 4c) are 
plotted in Fig. 9. In two areas to the left and right of the web (0.0 £ 0.4 and 0.7 ^ s ^ 

1 .0), the results are in good agreement. Locally, near the web (0.4 £ sz 0.7) the 
distributions differ. The peak values computed for the local solid model that included the 
web and noodle (black triangles) exceed the reference values (blue circles) by about 51%. 
The peak location is also offset. 



The results computed for a larger transition radius (r=5.08 mm as shown in Fig. 4c) are 
plotted in Fig. 10. In the area to the left of the web (0.0 s 0.4), the results are in good 

agreement. For the local solid model that included the web and noodle (black triangles), an 
additional local maximum was observed for s ~ 0.4 (G r / G c 5) which was not observed in 
the results used as reference (blue circles). In the areas right of the web (0.7 s s s 1.0), the 
results are higher for the models which included the web and the noodle. For this transition 
radius, the peak value is almost identical to the peak failure index computed for the 
reference results. 

For closer comparison, the computed failure index distributions across the width of the 
stringer are plotted in Fig. 1 1 and 12 for all models discussed. In the area to the left of the 
web (0.0 £ ££ 0.4), the failure index is low and all results are generally in good agreement. 
The peak values computed for the local solid models that included the web and the models 
that included the web and the noodle exceed the reference values (blue circles), which 
were obtained from models where the web had been modeled with shell elements. The 
location where the peak failure index was observed shifted from model to model. The 
results obtained from models where only the web was included in the local solid model 
(red squares) were in excellent agreement with the failure indices computed from the 
model with the smallest radius (r=0.254 mm, black triangles) and, therefore the smallest 
noodle region (Fig. 1 1). In two areas to the left and right of the web (0.0 £ 0.4 and 0.7 £ 

s £ 1.0), the results are in good agreement with the reference results. Locally, near the web 
(0.4 £ ££ 0.7), the distributions differ up to 67% as discussed above. For models with 
larger radii (r=0.71 1 mm, green diamonds; r=2.54 mm, orange triangles and r=5.08 mm, 
red squares), an additional local maxima or plateau was observed for 0.4 < s < 0.5 as 
shown in Fig. 12. With increasing radius, the peak value decreases and shifts to the right of 
the center 0.5 < s < 0.6. Also, in the area to the right of the web (0.7 £^^1 .0), the results 
are lower for the models with larger transition radius. 



5 Summary and Concluding Remarks 

The skin/stringer separation of a graphite/epoxy composite panel reinforced with three 
stringers and subjected to shear loading was studied using computational fracture 
mechanics analysis. The shear loading causes the panel to buckle and the resulting out-of- 
plane deformation initiates skin/stringer separation at the location of an embedded defect. 
The panel and surrounding load fixture were modeled with shell elements. 

A small section of the stringer foot, web and noodle as well as the panel skin in the vicinity 
of the delamination front were modeled locally with 3D solid elements. The mixed-mode 
strain energy release rates were calculated along a straight delamination front across the 
width of the stringer foot using the virtual crack closure technique. A failure index was 
calculated by correlating the results with a mixed-mode failure criterion of the 
graphite/epoxy material. Computed failure indices were compared to corresponding results 
where the entire web was modeled with shell elements and only a small section of the 
stringer foot and panel was modeled locally with solid elements. 

The results showed the following 

• including the stiffener web in the local 3D model increased the computed failure 
index by about 49% compared to the reference results where the web was modeled 
with shells 

• including the web, the noodle and transition radius in the local 3D solid model 
increased the magnitude of the failure index up to 67% and changed the local 
distribution across the width 

• the magnitude of the failure index decreased with increasing transition radius. 
Differences in the failure index distributions were likely caused by a different local 
deformation behavior due to different local stiffnesses of the models studied. The local 
differences in stiffness were mainly caused by the local modeling of the noodle and the 
transition radius. Based on the increase in computed failure index, it is suggested to use a 
high fidelity model including the noodle and transition radius whenever accurate analysis 
results are required. The results of this study are intended to be used as a guide for 
conducting finite element analyses of structures such as stiffened panels. In particular, this 



guidance is aimed towards analyses that attempt to simulate delamination growth and 
debonding. 
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Figure Captions 

Figure 1. Buckling of stringer stiffened panel subjected to shear loading (1016 mm x 
1016 mm) 

Figure 2. Local 3D model of panel skin and stiffener foot 

(a) Model with refined mesh near the edges (BF-8/10/8) 

(b) Model with refined mesh near the center (BF-20/20/20) 

Figure 3. Local 3D model of panel skin, stiffener foot and web 

(a) Model with refined mesh near the center (TFL3 -20/20/20) 

(b) Model with refined mesh near the edges (TFL3-8/10/8) 

Figure 4. Local 3D model of panel skin, stiffener foot and web including the noodle and 
transition radius 

(a) Model with refined mesh near the edges (TN1L3C-8/10/8) 

(b) Model with refined mesh near the center (TN1L3C-20/20/20) 

(c) Detail of radius and noodle for models with different radii 
Figure 5. Mixed-mode failure criterion 

Figure 6. Computed failure index - 3D web results versus reference solution 
Figure 7. Computed failure index - results for transition radius r=0.71 1 mm versus 
reference solution 

Figure 8. Computed failure index - results for transition radius r=2.54 mm versus 
reference solution 

Figure 9. Computed failure index - results for transition radius r=0.254 mm versus 
reference solution 

Figure 10. Computed failure index - results for radius r=5.08 mm versus reference 
solution 

Figure 11. Comparison of results for 3D web, radius r=0 .254mm and reference solution 
Figure 12. Comparison of results for different radii and reference solution 
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Figure 1 . Buckling of stringer stiffened panel subjected to shear loading (1016 mm x 1016 mm) 




(a) Model with refined mesh near the edges (BF-8/10/8) 



(b) Model with refined mesh near the center (BF-20/20/20) 
Figure 2. Local 3D model of panel skin and stiffener foot 




(b) Model with refined mesh near the edges (TFL3-8/10/8) 
Figure 3. Local 3D model of panel skin, stiffener foot and web 
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(a) Model with refined mesh near the edges (TN1L3C-8/10/8) 


radii r= 0.254 mm, 0.711mm, 
2.54mm, and 5.08mm modeled J 
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(b) Model with refined mesh near the center (TN1L3C-20/20/20) 



Figure 4. Local 3D model of panel skin, stiffener foot and web including the noodle and 
transition radius 
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Figure 5. Mixed-mode failure criterion 
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Figure 6. Computed failure index - 3D web results versus reference solution 
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Figure 7. Computed failure index - results for transition radius r=0.71 1 mm versus reference solution 
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Figure 8. Computed failure index - results for transition radius ^2.54 mm versus reference solution 
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Figure 9. Computed failure index - results for transition radius 7=0.254 mm versus reference solution 
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Figure 10. Computed failure index - results for radius 7=5.08 mm versus reference solution 
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Figure 1 1 . Comparison of results for 3D web, radius r=0.254mm and reference solution 
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Figure 12. Comparison of results for different radii and reference solution 



